Impurity pinning in transport through ID Mott-Hubbard and spin gap insulators. 
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A low energy crossover induced by Fermi liquid reservoirs in transport through a ID Mott- 
Hubbard insulator of finite length L is examined in the presence of impurity pinning. Under the 
assumption that the Hubbard gap 2M is large enough: M > Tl = Vc/ L (vc- charge velocity in 
the wire) and the impurity backscattering rate Fi <^ Tl, the conductance vs. voltage/temperature 
displays a zero-energy resonance. Transport through a spin gapped ID system is also described 
availing of duality between the backscattered current of this system and the direct current of the 
Mott-Hubbard insulator. 
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Recent developments in the nano-fabrication technique 
have allowed to produce relatively clean one channel 
wires of 1 — IQjjLm length. In these wires, ID Tomonaga- 
Luttinger Liquid (TLL) behavior |^ was observed in 
transport measurements |^,^. Then it has been sug- 
gested that the correlated insulating behavior may be 
tuned with this setup This behavior of the ID elec- 
tron systems is expected Q at half-filling where any 
Umklapp scattering has to open a Hubbard gap 2M in 
the charge mode spectrum of the infinite wire if the for- 
ward scattering is repulsive. Earlier consideration pre- 
dicted crossover from Fermi liquid transport carried by 
the reservoirs electrons at low energy to that of a Mott- 
Hubbard insulator in the one channel wire of long enough 
length L: M > Tl = Vc/L {vc- charge velocity in the 
wire). Charge in this insulator is carried by soliton exci- 
tations of the condensate emerging in the charge mode. 
This insulator, which is probably the easiest for the ex- 
perimental observation, features two marginal quantities 
0: charge of the soliton is unchanged e {e = fi = 1 
below) and the exponent 1/2 characterizing the electron- 
soliton transition brings about similarity with the free 
electron tunneling through a resonant level. 

In this paper we investigate the effect of impurity 
backscattering of low rate Fi <C M on the above 
crossover. For low energy following Schmid Q the prob- 
lem is mapped through a Duality Transform onto the 
model of a point scatterer with pseudospin imbedded in 
TLL. It is solved exactly by fermionization. The re- 
sult shows that the impurity enlarges the width of a 
zero energy resonance in the conductance up to Fi -I- F2 
where the exponentially small F2 (F2 ~ VT^Me"^*^/^^ 



at half filling) is the rate of tunneling of the condensate 
phase. The linear bias conductance at zero temperature 



G = 



reveals an exponential enhancement of the 



amplitude of the impurity potential yTi/T^ by the con- 
densate tunneling. This suppression of the resonant con- 
ductance does not affect the saturation of the current at 
J = F2 above the crossover as far as Fi -I- F2 ^ Ti,. Fi- 
nally, the result is addressed to the case when the gap 
opens in the spin mode of the wire. It occurs if the ef- 
fective interaction between electrons inside the wire is 
attractive Although it is not expected for the semi- 
conductor nanostructure, it has been observed in quasi 
ID organics ||]. The current backscattering in this sys- 
tem coincides with the direct one of the Mott-Hubbard 
insulator multiplied by a factor ^ , where F^ has the same 
meaning for this system as F; above. So, the backscatter- 
ing current is zero in the absence of impurities, and the 
current is maximum V/t:. Impurity appearance gives rise 
to backscattering together with simultaneous tunneling 
between the spin mode vacua of the wire. The conduc- 
tance is suppressed below a crossover energy and recovers 
above it. It results in a constant shift AJ = —r[ of the 
current at high voltage. 

Transport through a ID channel wire confined be- 
tween two leads could be modeled by a ID system 
of electrons whose pairwise interaction is local and 
switched off outside the finite length of the wire. Apply- 
ing bosonization and spin-charge separation we can de- 
scribe the charge and spin density fluctuations Pbix, t) = 
{dx<t)b{x,t)) / b ^ c,s, respectively, with (charge 
and spin) bosonic fields (j)c,s- Without impurities their 
Lagrangian symmetrical under the spin rotation reads 
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where ip{x) — 9{x)6{L — x) specifies a one channel ity(energy) in the channel. The parameter /ic = /x varies 
wire of the length L adiabatically attached to the leads the chemical potential inside the wire from its zero value 
X > L,x < 0, and vf{Ef) denotes the Fermi veloc- at half-filling and = 0. The constants of the forward 
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scattering differ inside tlie wire gb{x) — gt for x G [0,L] 
from those in the leads gb{x) = 1, and an Umklapp 
scattering (backscattering) of the strength Uc{Us) is in- 
troduced inside the wire. The velocities Vc^six) change 
from vp outside the wire to some constants Uc,s inside 
it. We can eliminate them rescaling the spacial coordi- 
nate Xoid in the charge and spin Lagrangians of (|lj) into 
Xnew = /q dy/vc,s{y)- As a result, the new coordinate 
will have an inverse energy dimension and the length of 
the wire becomes different for the charge mode L — > 1/Tl 
and spin mode L —^ Applying renormalization- 

group results of the uniform sin-Gordon model at ener- 
gies larger than T]^ or we come to the renormalized 
values of the parameters in (|l|). For the repulsive in- 
teraction when initially gs > 1 > gc-, the constant Us of 
backscattering flows to zero and 5's to 1, bringing the spin 
mode into the regime of the free TLL . The constant Uc 
of Umklapp process increases, reaching vp/vc at the en- 
ergy cut-off corresponding to the mass of the soliton M 
if the chemical potential /x is less than M. Meanwhile, gc 
flows to its free fermion value gc = 1/2 . These values 
of the parameters specifies the Mott-Hubbard insulator. 
For the attractive interaction initially gs < 1 < gc, and 
both constants C/c,s flow vice versa resulting in the spin 
gap insulator and TLL of some > 1 in the charge 
mode. Both cases of the interaction remain to some ex- 
tent symmetrically conjugated under the spin-charge ex- 
change even after accounting for a weak backscattering 
on a point impurity potential inside the wire < xq < L: 

2Vimp (t)c{t,Xc) (t>s{t,Xs) , , 

Cimp = ^^cos( p h(^o)cos( p ) (2) 

7ra V2 V2 

where Xc.s — xq/vc^s, = y^+fJ'Xc includes a phase of the 
scatterer (p. The amplitude of the potential Vimp spec- 
ifies transmittance coefficient [|o| as 1/(1 -I- V^^p), and 
a ~ i/Ep is momentum cut-off assumed to be deter- 
mined by the Fermi energy. Below we will first elaborate 
the case of the Mott-Hubbard insulator adjusting results 



Here 5o[0s] = dr ^ dx{{dr4's{T,x)Y + 

{dx4's{TTx)Y} / {%tt) is the free TLL Eucleadian ac- 
tion. All r-integrals run from to inverse temperature 
(3 = 1/T and aj = 0. To have all iTi^a-matrices time- 
ordered under the sign T, we attributed each of them to 
a corresponding time t assuming that their time evolu- 



to the transport through the spin gap insulator at the 
end. 

Duality Transform - An effective model for energies 
lower than some cut-off D' specified below may be read 
off following Schmid |8| from the expression for the par- 
tition function Z associated to the Lagrangian (|l]) plus 
(^). Without impurities the spin and charge modes are 
decoupled. After integrating out (j)c in the reservoirs the 
charge mode contribution into Z describes rare tunneling 
between neighbor degenerate vacua of the massive charge 
mode in the wire characterized by the quantized values 
of \/2'/>c (t, x) + 2 fix = 27rm, m is integer. Variation of m 
by a = ±1 relates to passage of a (anti)soliton through 
the wire ((anti)-instanton in imaginary time t). The tun- 
neling amplitude has been found |l[] by mapping onto a 
free fcrmionic model to be Pe"""/"'"^ , sq = y^M^~^Ji? = 
M sin -TO, [fi <C M) with the pre-factor P and the energy 
cut-off D' proportional to on approaching the pertur- 
bative regime ^ M. The calculation |0] by instanton 
techniques has corr ected P = C X \/i5'(sin^ roMTt)!/'* 
with the constant C of the order of 1. The parameter 
D' is a high-energy cut-off to the long-time asymptotic s 
of the kink-kink interaction: F(r) = IhI^/t^ + l/D'^} 
created by the reservoirs. It varies with /i from D' ~ 
y/MTE at /i to £>' ~ {M/h)Tl ff > Tl- A 
crucial modification to this consideration produced by 
the impurity under the assumption EpVimp <C M en- 
sues from the shift of the m-vacuum. Since it is equal 
to (— 1)™^^^ cos(/3COs(^^^^j^) the neighbor vacua be- 
come non-degenerate. This can be accounted for by as- 
cribing opposite values of the pseudospin variable cr = ±1 
to the neighbor vacua which are the eigenvalues of the 
third component of the Pauli matrix (73. The energy 
splitting becomes an operator (73^^^ cos (/s cos (^^^^i^) 
acting on the pseudospins, and every (anti-)instanton 
tunneling rotates a a^-vahie into its opposite with the 
Pauli matrix tri. The partial function then can be writ- 
ten as 



(3) 



tion is trivial. Noticing that the interaction F coincides 
with the pair correlator of some bosonic field 9c whose 
evolution is described with the free TLL action 5o[6'c], 
we re-write (||) ascribing a factor exp{^9c{Tj, 0)/\/2) to 
the Tj (anti-)instanton, respectively: 
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It is easy to recognize a standard Hamiltonian form Z, — est x Tr{e ^'^} in with 

H = 7io[<^s(a;)] + 7^o[ec(^)] - 2Pe-'°/^-a, cos(0,(O)/V2) - cos(^)a3 cos(^%l) (5) 



Here (l)s{x) and 0c (a;) are Schrodinger's bosonic op- 
erators related to the variables 0s (j, a;) and 6c{t,x) 
of the functional integration in fl). The opera- 
tor Ho[(j)s{x)] {Hq{6c{x),Vc)), a function of the field 
(j)s{x) {Oc{x)) and its conjugated is a free TLL Hamil- 
tonian (g = 1) corresponding to the free TLL action 
5o[(/)s] {Sq\9c]) in (Q), respectively. The Dual model spec- 
ified by (||) is equivalent to the initial one (||) at low en- 
ergy. It relates to a Point Scatterer with internal degree 
of freedom in TLL. Fortunately, this in general rather 
complicated model may be solved easily through fermion- 
ization in our particular case. This simplification stems 
from the marginal behavior of the Mott- Hubbard insula- 
tor: the charge of the transport carriers does not change 
on passing from the low energies to the higher ones de- 
spite the nature of the carriers does. 

Fermionization - From the commutation relations and 
hermiticity, the Pauli matrices can be written as CTq — 
e"'^'''?/3C7 witli Majorana fermions ^1,2,3 
and antisymmetrical tensor e : e^^^ = 1. Since the inter- 
action in (j^) is point-like localized and its evolution in- 
volves only the appropriate time-dependent correlators, 



we can fermionize it making use of: 

, , fly ,- g<=(o) 

MO) = \ TT^i^ ^ ■ 

V ZTT 

Here ipc,s{0) is the a; = value of the charge (spin) 
fermionic field, respectively. These fields have linear dis- 
persions taken after the related bosonic fields with mo- 
mentum cut-offs (equal to the energy ones) D' and Ep, 
respectively. Substitution of these fields into (||) produces 
a free-electron Hamiltonian where the interaction reduces 
to tunneling between the ipc.s fermions and the Majorana 
one ^2 (= below). Application of a voltage V between 
the left and right reservoirs forces us to use the real-time 
representation. Since each instanton tunneling transfers 
charge A0c/(\/27r) — le between the reservoirs the volt- 
age may be accounted for by a shift 6c/ \/2 -—^ 9c/ \pl -f Vt 
of the cos-argument in the real-time form of the action 
from Eq.(^. Assuming that it remains small enough, 
y < Ti < Af , we will neglect its effect on the other pa- 
rameters. The real-time Lagrangian associated with the 
fermionized Hamiltonian (^ reads: 



+ h.c] - Vri[Mio,mt) + h.c] 



(6) 



where the rate of impurity scattering is Fi — 
{cos{tp)Vi„ip)^ and the rate of the instanton tun- 

nehng is T2 = 2ttC'^ VTlM sin^ we''^^^/'^'^ . The cur- 
rent flowing through the channel is J = ~^f% — 

-iv^[^/'+(0,t)^(i)e~'^* - h.c.]. Its calculation with 
the non-equilibrium Lagrangian ^ needs avail of the 
Keldysh technique. Transforming voltage into the non- 
zero chemical potential of the fermions equal to —V, we 
can write in standard notations 



J = 2WT 



^ReG>Ju;) 
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T 



Gf{u)+G>{u) 



(7) 



where / is the Fermi distribution function. To find the 
Green functions G^, of the Majorana field, its free 

2/(cj±i0), = ±2m5{uj) 



Green functions :(7^*'^'' — 



are substituted into the appropriate Dyson equations 

qR(A) ^ gR(A){l + Y.^{A)qR{A)^ 

G^ = (1 + GV)^^! + S-^G^) + G^E''G^ (8) 

where 7 stands for <,>,c, c. The self-energies for the 
Majorana field are:E^^'^^ = ipi(ri-|-r2), = -i(r2(2- 
fiiuJ ± y)/T)) + 2(1 - fiuj/T)). Theh substitution 
into (||) and into (|^), subsequently, results in the current: 



J = 



2F2(ri + r2) f /(^)-/(^) 



(9) 



which is the current passing through a resonant level of 
the half-width 2(Fi + T2) and suppressed by the factor 
r2/(ri + F2). The typical features of this current can be 
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illustrated with its zero temperature behavior vs. volt- 
age and the linear bias conductance G vs. temperature, 
respectively: 



2V. 



■ arctan 



V 



2(ri- 
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(10) 
(11) 



where ip' {x) is the derivative of the di-gamma function, 
'0'(l/2) = 7r^/2, and the high temperature asymptotics 
of G is r2/(2T). The zero-temperature conductance 
G ^(rf+r2) i^a-nifest in (^) is the function of Fi /T2 (x 

(cos (pVimpe'^''^'^^ )'^Ef I \/ TlM siv? w. Comparison with 
the initial transmittance of the impurity scatterer shows 
that ^jTi/T2 may be conceived as a renormalization of 
the initial amplitude Vimp by the instanton exponent and 
a power factor of the TLL for g = 1/2. 

Spin gap insulator - In this case the chemical poten- 
tial of the spin mode always lies in the center of a gap 
2M. Meanwhile the charge mode is inhomogeneous TLL 
with the constant gc = g > l((7c = 1) inside (outside) 
the wire. Under Duality Transform applied to the mas- 
sive spin mode the partition function Z' takes the form 
of Eq.(^ with the following modifications. The field 
(j)s {9c) is changed by 4>c {Os), respectively, and Xg by 
Xc- The factor cosip disappears, as is accumulated by 



0c- The amplitude of the instanton tunneling between 
the spin mode vacua becomes P'e^*^/^^-, where a new 
pre-factor P' and cut-off D" are specified by the same 
expressions which were written for P and D' , respec- 
tively, after substitution of instead of into them. 
The action for the (pc field is inhomogeneous: S[(j)c] — 
Jl^dT J dx{{drMT,x))^ + id,Mr,x)y}/{8gc{x)n). To 
proceed in line with the above scheme of calculation, we 
have to first integrate out the quick modes of the (jjc 
field whose energies are larger than T^. Unless Vc/vs 
becomes essentially small at g ^ 1, these energies lie 
much higher than the expected crossover. Therefore, 
we can integrate neglecting the weak tunneling between 
the spin mode vacua. In the lowest order in vari- 
ation of the cut-off from Ep to results in multipli- 
cation of V^^p by a isiCtov:{EF /TLfs F[xcTl) where 
F{xcTl) accounts for the interference produced by the 
inhomogeneity of the forward scattering. This func- 
tion may be extracted from the lon g-ti me asymptotics 
< ^i^^(x,t)^-i<i>,(xfi) > calculated |[l|] with 5[(^c] as: 

F2(z) =COnst(z2+^2)<,r-Q-^^J(^±^)2^^2]gr^™±i ^ 

consi'x [(z2+72)((l-z)2-f72)]9'-, where 7 = Tl/Ef and 
r = ijp^ is a negative reflection coefficient. At the ener- 
gies below the new cut-off Tl , the inhomogeneous action 
S[4>c] uiay be changed by the homogeneous one iSo[^ic]- 
Then exercising the above fermionization we come to a 
real-time Lagrangian analogous with (§: 



i^dt^ + i 



dxi^t{dt + d,)i>a - Vr'i[^+(0,t)e(i)e-*''' + h.c] - Vr^[^,+ (0, t)C(i) -f h.c 



(12) 



Here the rate of the instanton tunneling T'2 = 
y'T^Me"^*^/-^^ has an manifest similarity with 
r2 and the rate of impurity scattering T'l = 
^{EF/TLy-9F{xcTL)V^^p acquires an additional 

smaU factor {Ef /TlY^^ F{xcTl), as gc > 1. Creation 
and annihilation of the fermion of the charge mode in the 
fermionized model describes a change of chirality of one 
quasiparticle inside the wire. Therefore, the term pro- 
portional to in ( p^ leads to the backscattering cur- 

rent: J^sc = = -ty/T[[t^+ {0,t)at)e-^''' - h.c] 

equal to minus deviation of the direct current from its 
maximum V/n value. The new Lagrangian ( p^ and the 
current Jbsc transform into the old Lagrangian (^) and 
the current J, respectively, on changing F'^, r2 by Fi. 
Hence the average value of Jbsc may be extracted from 
(|) in the same way. In particular, the linear bias con- 
ductance: G = (1 - ^V' (1/2 + {F'l + F'2}/(7rT)))/7r is 
strongly suppressed at T = as G = ^p, ^ . Similarly 
to the Mott-Hubbard case it shows that the initial am- 
plitude of the potential Vtmp renormalizes to \/^i/^'2 oc 

V.mpe^'/^^ ^ Ef{Ef/TlY-S'F{xc)I ^/T[M at low en- 



ergy. Above the crossover temperature 2{T'i -I- Fj) the 
conductance recovers as ^ ^1 — . The backscatter- 
ing current related to ( |lO| ) results in the constant shift 
of the direct current J = V/n — T'l above the voltage 
crossover 2{T{ -I- F^. 

Finally, we have shown that a zero-energy resonance 
in transport through a ID Mott-Hubbard insulator of fi- 
nite length predicted [Q for the clean wire stands weak 
impurity pinning while the rate of the impurity backscat- 
tering Fi is small Fi <C T^. We have also described how 
opening a gap in the spin mode suppresses the charge 
transport. 
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